addpath('npy-matlab-master\npy-matlab')
A1 = readNPY('map_MCD_30K_niI2_test_x.npy')
A1=A1*1e5/2
surf(A1,'edgecolor','none')
set(gca,'xlim',[1 41],'ylim',[1 41])
view(2)
%colormap(winter)
caxis([-1.3,1.3])
set(gca,'box','on','Xcolor',[0 0 0],'Ycolor',[0 0 0])
xticks({})
yticks({})
colorbar('Ticks',[-1.2:0.6:1.2],'TickLabels',{})
set(gcf,'Position',[0, 0, 380, 300])